'Magic Melters' Have Geometrical Origin. 
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Recent experimental reports bring out extreme size sensitivity in the heat capacities of Gallium 
and Aluminum clusters. In the present work we report results of our extensive ah initio molecular 
dynamical simulations on Ga3o and Ga3i , the pair which has shown rather dramatic size sensitivity. 
We trace the origin of this size sensitive heat capacities to the relative order in their respective 
ground state geometries. Such an effect of nature of the ground state on the characteristics of heat 
capacities is also seen in case of small Gallium and Sodium clusters indicating that the observed 
size sensitivity is a generic feature of small clusters. 

PACS numbers: 31.15.Qg,36.40.Sx,36.40.Ei,82.20.Wt 



The finite temperature behavior of clusters has shown 
many interesting and intriguing properties 0, 0, 0, . 
Recently the calorimetric measurements reported by Jar- 
rold and coworkers found that small clusters of Tin and 
Gallium in the size range of 17-55 atoms have higher than 
bulk melting temperatures (Tni[ bulk ]) 00. A striking ex- 
perimental result from the same group showed extreme 
size sensitivity in the nature of the heat capacity for Ga 
clusters in the size range of 30-55 atoms 4J. It turns 
out that addition of even one atom changes the heat ca- 
pacity dramatically. For example, Ga^ has a rather flat 
specific heat curve whereas heat capacity of Ga^ has a 
well defined peak and has been termed as 'magic melter'. 
A similar size sensitive feature has also been observed in 
the case of Al clusters p. 

The explanation and understanding of various experi- 
mental observations have come from the first principles 
Density Functional (DF) simulations 0,0,0, 110- For 
example, the higher than bulk melting temperature for 
Sn and Ga clusters is understood as due to the differ- 
ence in the nature of bonding between the cluster and 
the bulk 0,0,0- However, the extreme size sensitiv- 
ity displayed in Gallium and Aluminum clusters is still 
an unexplained phenomena. The present work addresses 
this issue by employing first principles DF methods. In 
this letter we report our results of ab initio molecular 
dynamical (MD) simulations carried out on Ga3o and 
Ga3i . It is of some interest to note that similar size sen- 
sitive heat capacities have been observed in case of Ga n 
(n = 17, 20) and Na n (n = 40, 50, 55) [ll| clusters. In 
both these cases addition of few atoms changes the nature 
of heat capacities significantly. By analyzing the geom- 
etry of the ground state, we establish a definitive corre- 
lation between the nature of the ground state and the 
observed heat capacity. Our detailed calculations show 
that an 'ordered' ground state leads to a heat capacity 
with a well defined peak while a cluster with 'disordered' 
ground state leads to a flat heat capacity with no dis- 
tinct melting transition. In what follows we will make 
the meaning of 'order' and 'disorder' precise and provide 
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FIG. 1: The heat capacity of Ga3i, Ga3o computed over 
90 ps. 



explanation for the size sensitive heat capacities. 

We have carried out isokinetic Born-Oppenheimer MD 
simulations using ultrasoft pseudopotentials within the 
generalized gradient approximation (GGA) For all 

the clusters reported here we have obtained at least 200 
equilibrium structures. For computing heat capacities 
of Ga3o and Ga3i the MD calculations were carried out 
for 16 different temperatures, each with the duration of 
150 ps or more, in the range of 100 <T< 1100 K, which 
results in a total simulation time of 2.4 ns. In order 
to get converged heat capacity curves especially in the 
region of coexistence, more temperatures were required 
with longer simulation times. The resulting trajectory 
data has been used to compute the ionic specific hea t by 
employing the Multiple Histogram (MH) method [lllll. 

Fig. shows the calculated heat capacity of Ga3o and 
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FIG. 2: The ground state geometry of Ga3o and Ga3i with 
two different perspectives, perspective (b) is rotated by 90° 
with respect to perspective (a). 



Ga3i. Evidently the dramatic difference in the heat ca- 
pacities of Ga3o and Ga3i, observed in the experiments 
is well reproduced in our simulations. Thus Ga3i has 
a well defined peak in the heat capacity whereas the 
heat capacity for Ga3o is rather flat. We also note that 
both Ga3o and Ga3i becomes 'liquid-like' at tempera- 
tures much higher than Tm[ bulk ] (303 K) i.e. around 
500 K, consistent with the experiments. In order to gain 
insight in to these observations we analyze the ground 
state of Ga3o and Ga3i. In Fig. [21 we show the ground 
state geometries of Ga3o and Ga3i with two different 
perspectives. A cursory analysis of Fig. 0-a may lead 
to the conclusion that the only difference between Ga3o 
and Ga3i ground states is the presence of the capped 
atom in Ga3i. However, a different view obtained by 
rotating the cluster by 90° brings out the significant dif- 
ferences in Ga3o and Ga3i, clearly indicating that Ga3i 
is more ordered. A careful examination of Fig. El b shows 
the presence of well ordered planes in Ga3i. Such planes 
are only in a formative stage and considerably deformed 
in Ga3o- In fact an addition of just one atom in Ga3o 
displaces all the atoms by a significant amount which 
makes Ga3i more ordered. That a single atom makes a 
substantial change is also seen by the fact that there is 
noticeable difference in the coordination number in these 
two clusters. In Ga3o, 5 atoms have 4 or more coordi- 



nation number whereas in Ga3i, 14 atoms have 4 fold 
or more coordination. Therefore we termed Ga3o as a 
'disordered' structure relative to Ga3i. 

Thus when the system is disordered or amorphous each 
atom (possibly a group of atoms) is likely to have differ- 
ent local environment. That means different atoms are 
bonded with the rest of the system with varying strength. 
Consequently, their dynamical behavior as a response to 
temperature will differ. Some of the atoms may pickup 
kinetic energy at low temperatures while the others may 
do so at higher temperatures. In a given structure if a 
large group of atoms are bonded together with a similar 
strength forming an island of local order it is reasonable 
to expect that they will 'melt' together. In this case the 
cluster can be considered as (at least partially) ordered 
and will show a well defined peak in the heat capacity. 
However, if the system is disordered in the sense that 
there are no such islands of significant sizes having lo- 
cal order then we expect a very broad continuous phase 
transformation. Indeed, our analysis of mean square dis- 
placements (MSDs) for individual atoms brings out this 
fact clearly. The MSDs for individual atoms is defined as 

1 M 

< r ? W> = mE [ r '(W + t)- R/(W)] 2 , (1) 

m=l 

where R/ is the position of the Jth atom and we aver- 
age over M different time origins to m spanning the entire 
trajectory. In Fig.[2|we show MSDs of individual atoms 
for Ga3o and Ga3i at 250 K. The contrast between the 
kinetic response of individual atoms in Ga3o and Ga3i is 
very clear. For Ga3o, the MSDs of individual atoms show 
that some of the atoms (at least 10) have picked up more 
kinetic energy compared to others and hence have signif- 
icantly higher displacements (9.0 A 2 as compared to 0.45 
A 2 ) whereas in Ga3i all atoms are oscillating about their 
mean positions and exhibit small values of MSDs (0.45 
A 2 ). Thus MSDs clearly indicate that in Ga3o different 
atoms have different mobilities. This wide distribution 
of MSDs in Ga3o indicates that the cluster is in coexis- 
tence phase around 250 K and is continuously evolving. 
This is precisely what is expected if the cluster is disor- 
dered in the sense described above. This phenomena has 
also been observed in the extended systems and a simi- 
lar analysis has been used to characterize the nature of 
spatial inhomogeneities with considerable success, yjj 

The difference in the mobilities of individual atoms 
in these two clusters is also reflected in the root-mean 
square bondlength fluctuations (5 rms ) shown in Fig. 
S rms shows a clear signal for the beginning of the change 
of phase around 450 K for Ga3i. However, in case of 
Ga3o the transition is spread over a much broader range 
of temperatures. In fact the coexistence region for Ga3i 
is over 175 K (from 425 K to 600 K) and for Ga 30 it 
extends over 425 K. It is interesting to note that in Ga3o 
the isomerization begins around 175 K and continues till 
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FIG. 3: The MSDs for individual atoms of Ga3o and Ga3i 
computed over 90 ps. 
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FIG. 4: The S r ms for Ga3o and Ga3i computed over last 90 ps. 




FIG. 5: The Electron Localization Function of (a) Ga3o and 
(b) Ga3i at the value of xelf = 0.68. The black lines show 
connected basins. 



600 K. 

The nature of the 'order' can also be brought out by ex- 
amining the Electron Localization Function (ELF). ELF 
has been found to be extremely useful for elucidating the 
bonding characteristics of a variety of systems . For a 
single determinantal wavefunction built from KS orbitals 
ipi, the ELF is defined as, 



where 



XELF = [1 + {D/D h f\ 



D h = (3/10)(37r 2 ) 5/ V/ 3 , 



(2) 



(3) 



D = (1/2) £ |V^| 2 -(1/8)|V P | 2 /P, (4) 

i 

with p = p(r) is the valence-electron density. ELF is de- 



fined such that values of Xelf approaching unity indicate 
a strong localization of the valence electrons and covalent 
bonding. ELF partitions the molecular space into regions 
or basins of localized electron pairs or attractors. At very 
low values of the ELF all the basins are connected. In 
other words, there is a single basin containing all the 
atoms. As the value of the ELF is increased, the basins 
begin to split and finally we will have as many basins as 
the number of atoms. The value of the ELF at which the 
basins split is a measure of the strength of interaction 
between the different atoms. Fig. shows the isosurface 
of ELF taken at Xelf = 0.68, for Ga3o and Ga3i. It can 
be noted that for Ga3i, 26 atoms are connected via a 
single basin whereas for Ga3o the largest basin contains 
12 atoms with other 'fragmented' basins. This supports 
our earlier observation that Ga3i has significantly more 
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FIG. 6: The heat capacity of Na4o, Naso and Nass computed 
over 90ps data. 

similarly bonded atoms than Ga3o- Further evidence for 
the disordered or amorphous nature of Ga3o comes from 
the comparison of entropies of these systems (Fig. not 
shown). As expected the entropy of amorphous structure 
(Ga3o) rises rather sharply as compared to Ga3i (which 
is more ordered). Quite clearly, the amorphous nature 
leads to substantially large number of accessible states in 
case of Ga3o and is more by a factor of ten as compared 
to Ga3i in the low energy region. 

As mentioned earlier this size sensitive behavior is not 
unique to the Gallium clusters reported here and has 
been observed in small Ga 9], Al p, and Na [ll| clus- 
ters. As an example we show the heat capacities of Na n 
(n = 40, 50, 55) clusters in Fig.0 The change in the na- 
ture of heat capacity as the cluster grows from Na4o to 
Na55 is quiet evident from Fig. Our detailed analysis 
of the ground state geometries shows a direct correlation 
between the nature of the ground state and calculated 
heat capacities. It may be noted that Nass is highly 
symmetric and very well ordered. Na4o is also ordered 
and has basin containing substantially large number of 
atoms but Naso is relatively disordered which is clearly 
reflected in their heat capacities flT^ . 

The main contribution of the present work is to bring 
out a definitive relationship between the local order in the 
cluster and its finite temperature behavior. As the cluster 



grows in size it is very likely that it will evolve through a 
succession of such ordered and disordered geometries. In 
such cases addition of one or few atoms is likely to change 
(as demonstrated in this work) the nature of the ground 
state abruptly. Thus, the size sensitive nature of heat 
capacities is generic to small clusters and related to the 
evolutionary pattern seen in their ground states. The 
evidence for this comes not only from Gallium clusters 
but also from clusters of Sodium and Aluminum having 
very different nature of bonding. 
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